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We investigate the coupled dynamics of the local hydrodynamics and global electric response of 
an electrodialysis system, which consists of an electrolyte solution adjacent to a charge selective 
membrane under electric forcing. Under a DC electric current, counterions transport through the 
charged membrane while the passage of co-ions is restricted, thereby developing ion concentration 
polarization (ICP) or gradients. At sufficiently large currents, simultaneous measurements of voltage 
drop and flow field reveal several distinct dynamic regimes. Initially, the electrodialysis system 
displays a steady Ohmic voltage difference (AVohm), followed by a constant voltage jump (Af4). 
Immediately after this voltage increase, micro-vortices set in and grow both in size and speed 
with time. After this growth, the resultant voltage levels off around a fixed value. The average 
vortex size and speed stabilize as well, while the individual vortices become unsteady and dynamic. 
These quantitative results reveal that micro-vortices set in with an excess voltage drop (above 
AVohm + AUc) and sustain an approximately constant electrical conductivity, destroying the initial 
ICP with significantly low viscous dissipation. 


Ion concentration gradients emerge during a separa¬ 
tion process involving a charge selective surface (elec¬ 
trode or membrane) in an electrolyte solution, inducing 
a decreasing ion concentration towards the interface and 
thereby hampering ion transport. For example, for wa¬ 
ter purification using electrodialysis under electric forc¬ 
ing, the charge selectivity of an ion exchange membrane 
causes an enrichment of counter-ions on the permeated 
side of membrane and a depletion of co-ions on the feed 
side. This so-called ion concentration polarization (ICP), 
with a decreasing ion concentration on the (feed side) 
interface, is a common theme in electrochemical appli¬ 
cations that influence the performance of ion separation 
and transport. The diversity of ICP-associated applica¬ 
tions has recently motivated numerous studies, rangiM 
from micro-and-nano-junctions M , electrodialysis [J- 
[T^ , desalination [l3, El] , electrodeposition Eht - Elj , and 
fuel cells [^ . 

In electrodialysis, theoretical analyses of low dimen¬ 
sions reveal that ICP can drive hydrodynamic instability 
through an equilibrium or non-equilibrium (electro- 
osmotic or bulk electro-convective) mechanism I, [H- 
suggesting an additional charge transport due to 
ICP-induced fluid motion under a sufficiently large DC 
voltage. Recent advances have been made with direct 
numerical simulations (DNS) [28l - l^ . providing insights 
into ion concentrations and flow velocity adjacent to a 
charged membrane. Experimentally, under a pressure- 
driven micro-channel flow, the advection and height se¬ 
lection of the unidirectional sheared vortices were charac¬ 
terized along the membrane m , but the internal vortex 
structure was not probed. Furthermore, quantitative ex¬ 
periments of the flow field in electrodialysis without shear 
flow are still missing, in particular under high electrical 
forcing Q. In this paper, we show the first quantitative 


measurements of the coupled hydrodynamics and electri¬ 
cal response of an ionic solution in the vicinity of a charge 
selective membrane, under a constant DC electrical cur¬ 
rent without an external shear flow. 

Fig. [T] shows the experimental setup and the resulting 
fluid dynamics at high electric forcing in electrodialysis. 
We used a cation exchange membrane, CEM (Neosepta 
CMX, surface area of 3 mm x 4 mm with thickness of 
170 /rm), horizontally placed in a transparent PMMA 
cell filled with a lOmM CUSO 4 electrolyte between two 
copper electrodes. We performed chronopotentiometric 
measurements between the top cathode and the bottom 
anode with an electrometer (Autolab PCSTAT30 Poten- 
tiostat). This method consists of forcing a constant DC 
electric current, /, through the ionic solution across the 
membrane and measuring the time-dependent voltage 
difference, AU(t) between the electrodes. The bottom 
anode, where copper oxidizes, serves as a cation source; 
the top cathode, where copper reduces, acts as a cation 
sink. This configuration of the electrodes can suppress 
the occurrence of gravitational convection caused by a 
variation of fluid density due to ion-concentration since a 
heavier Cu^^-rich solution is present close to the anode 
at the bottom of the fluid cell [3l|. We noticed that Cu 
dendrites form on the cathode after long experimental 
runs (> 1000 s) and at a high current density. Thus, 
electrodes are cleaned before each experiment. We fo¬ 
cused on the time-series data before the maximal ICP 
condition for the cathode. 

In addition to the electrical measurements AV{t), si¬ 
multaneously, the flow motion is observed close to the 
membrane interface under different DC currents (see 
Fig. [T)D-d). The flow motions are observed by seeding 
0.1 wt% 5 p,m polystyrene tracer particles (Microparti¬ 
cles GmbH, with the particle density of 1.05 g/cm^) to 
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FIG. 1: (a) Schematic diagram of the experimental setup 
of a charge selective (cation exchange) membrane, which al¬ 
lows cations to pass the membrane while rejecting anions, im¬ 
mersed in an ion solution of CUSO 4 under an electric DC cur¬ 
rent between two copper electrodes. We simultaneously mea¬ 
sure the global voltage drop, AV, under a constant applied 
current, I, and local hydrodynamics of the ionic solution, (b- 
d) Representative motions of the seeding micro-particles re¬ 
veal different hydrodynamic states in time: (b) Electrical con¬ 
duction regime without hydrodynamic convection (at 100 s). 
This transforms to a convection regime with micro-vortices 
growing in size with time, (c-d) 220 s and 420 s, respectively. 
The applied DC current density, J, was 10 A/m^. The dashed 
lines indicate the bottom boundary of the membrane. 


the solution. For these nearly buoyancy-neutral micro¬ 
particles, the theoretical sedimentation speed is small 
(Ri 0.7 fimfs, estimated by the Stokes drag equation [s^l. 
The zeta potential of these micro-particles is measured 
to be « —1 mV for 1 mM CUSO 4 (Zetasizer Nano ZS, 
Malvern). These micro-particles do not influence the 
electrical response of the electrolyte solution, as shown 
by similar AV (t) data obtained with and without the 
micro-particles. To avoid particle aggregation, a non¬ 
ionic surfactant, TweenSO, (0.1 wt%) is added to the solu¬ 
tion. We obtain accurate flow fields employ ing a particle 
image velocimetry (PIV) technique (32|, |33|. The micro¬ 
particles are illuminated by a thin laser sheet (Firefly 
laser, 808 nm, Oxford Lasers) with a pulse duration of 
20 p,s and a pulse power of 0.3 mJ/pulse (illuminating 
3 mm X 2 mm x 200 fim). The scattered light is cap¬ 
tured at 20 Hz by a CCD camera (Sony XCG-H280E, 
1920x1080 px^), with a magnifying lens (Navitar, 2-14 
x) placed perpendicular to the laser sheet. Sets of 50- 
200 images are analyzed using Imaged software (NIH) to 
visualize and measure the vortex motions and sizes. 

We also determine the vortex speed and size using PIV 
analysis, with a typical time delay (of 0.1 sec) between 
the image pair. The focal depth of the optical system is 
measured to be « 200fim. Particles outside this depth of 
field are larger, have a lower light intensity, and are sys- 



FIG. 2: Chronopotentiometric measurements of the total elec¬ 
tric voltage drop, AV(t) changing in time across the cell, un¬ 
der different applied DG current densities, J, through the 
charge selective membrane. The total potential drop reveals 
three distinct responses: (1) The initial constant Ohmic re¬ 
sistance of the whole cell AVohm = I Rohm, (2) The critical 
voltage jump, A14 that starts at the transitional time, Tc, and 
(3) The over-limiting voltage drop over the vortex mixing re¬ 
gion, AVoic- 


tematically filtered in the PIV analysis. To calculate the 
flow field, we use a multi-grid cross correlation method 
with decreasing window size by w 50% [13, First we 
use a 128 x 128 px^ interrogation window to determine a 
reference vector field. This field is subsequently used to 
calculate a (window) shift for the next correlation. To get 
a higher resolution, the second calculation is done with 
windows of 96 x 96 px^, and the vectors are displayed 
with a 50% overlap [U. From these vector fields, we de¬ 
termine the mixing layer thickness, Lmix, of the vortex 
region and the root mean square velocity (average vortex 
speed) within this layer. 

Each experiment starts with a uniform concentration 
of CUSO 4 at both sides of the membrane. As a DC 
electric current is forced through the charged membrane, 
the counterions (Cu^"*" cations), which can easily pass 
through the membrane, migrate upwards, whereas co¬ 
ions (S 04 ^ anions) migrate downwards. Because the co¬ 
ions are retained by the charged membrane, the co-ion 
concentration enriches at the upper (cathode) side and 
depletes at the bottom (anode) side of the membrane. 
At a critical condition, a maximum ion-depletion occurs, 
with a vanishing co-ion concentration at the membrane 
surface. This so-called limiting current density can be 
estimated by balancing electromigration with diffusion 
of the co-ions [13, [111, and Jum = 3.3 Ajm^ for our 
electrodialysis system. We apply a current density, J, 
above this limit to study the ‘over-limiting’ conductivity 
induced by ICP, a long-standing unsolved problem. 

Fig.[2]shows the experimental results of the global elec¬ 
trical responses of the electrodialysis system. In time, 
the redox reactions at the electrodes start immediately, 
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and ion concentration polarization is gradually devel¬ 
oped. The initial constant voltages {AVohm) reflects a 
constant Ohmic resistance of the electrodialysis system: 
R = 11.2 fcO, obtained from a linear fit of AVohm for dif¬ 
ferent I. Subsequently, ion concentrations are depleted 
close to the anode side of the membrane. The depletion 
grows until the critical limiting condition, which leads to 
an increased electrical resistance. This critical condition 
is manifested by a sharp increase in measured voltage, 
AVc, after a transitional time Tc- By the same token for 
Jiim described above, can be estimated using Pick’s 
second law with a vanishing co-ion concentration (c ) at 
the membrane surface, i.e., Sand’s equation [l0|,[36|. In 
agreement with the theory, our Tc has a linear relation¬ 
ship with 1//^, with a fitted cation transport number in 
the membrane of 0.9, consistent with the previous ex¬ 
perimental results [13, m, [131 • The voltage jump AVc 
corresponds to the (electrical current) plateau region in 
the IV curve under a DC voltage (e.g. Fig. 4 in Ref. [13) 
and depends on the type of the membrane. We measured 
AVc = 0.3 ± 0.03 V for different J, ranging from 6.3 to 
16.7 A/m^. The jump is followed by a further voltage 
increase, AVoic, until AV reaches a saturated value at a 
later time. 

The first two characteristic electric responses, AVohm 
and the onset of AI 4 , are well understood, however the 
later-time AVoic and the transport mechanisms caus¬ 
ing over-limiting conductance, beyond AI 4 , have been 
extensively debated [l5j| . This phenomena is mani¬ 
fested in the additional voltage drop AVoic observed 
under a constant / or the increasing currents in the 
conventional current-voltage curves under DC voltages 
by other studies [H, HI]- Several mechanisms have 
been proposed, including water dissociation, hydrody¬ 
namic convection, and charge-induced membrane dis¬ 
charge [13, EE EE EE EE EE- Previous work has 

shown that water dissociation for our type of CMX mem¬ 
brane is insufficient to account for the observed over¬ 
limiting conductance EE El, EE- To gain insight, we 
analyze the coupled hydrodynamics from the captured 
images, using particle pathlines (Fig. [T] b-d) and PIV 
analysis (Fig. [3] b-d). 

In Fig. |3l a representative set of the coupled dynamics 
of electric response AV(t) and vertical vortex size Lmix 
(measured from the membrane surface) is shown. The er¬ 
ror bars represent the variation of the individual vortex 
sizes. Four distinct regimes are delineated in Fig. [3l The 
initial Regime I is electric Ohmic conductive, with no hy¬ 
drodynamic convection observed from the motion of the 
micro-particles. In Regime II, voltage jump AVc occurs, 
with an increasing electrical resistance of the electrodial¬ 
ysis cell. This starts at the transition time (e.g., Tc ~ 77 s 
for J = 10 A/m^). In Regime III, small counter-rotating 
vortex pairs appear along the membrane surface. The 
thickness of this mixing vortex layer Lmix grows linearly 
in time (e.g. Fig. [SJd). Simultaneously, AV{t) gradually 
grows in this regime. Finally, in Regime IV, both AV 
and vortex size saturate and fluctuate at fixed values. 


(a) 




FIG. 3: (a) Simultaneous measurements of the dynamics of 
vertical micro-vortex size (•) and voltage response AV (—), 
under a constant current density of 10 A/nV across the charge 
selective membrane. Four characteristic transport regimes are 
delineated by the dashed lines. I. electric Ohmic conduc¬ 
tion without hydrodynamic convection, II. a potential jump 
with the development of ICP, III. a hydrodynamic convec¬ 
tion regime with linearly growing electric resistance and vor¬ 
tex size in time, and IV. a saturated regime with saturated 
values of both vortex size and electric resistance, (b-d) The 
corresponding flow field, velocity vectors and vorticity (V x u) 
obtained with a PIV technique at different times; (b) and (c) 
show the growth and (d) the unsteady nature of the micro¬ 
vortices. The vertical arrow in (b-d) indicates the length 
scale of vortex mixing layer, Lmix, which initially increases 
with time. Lmix is measured to be 260, 480, and 550 /rm 
from (b) to (d), respectively, at different times {!) indicated. 

The vortices are observed to move laterally and merge to¬ 
gether, showing unsteady dynamics (see the supporting 
videos [IE). Consistent with our experimental findings, 
the unsteady and chaotic movements of saturated vor¬ 
texes have been observed in recent direct numerical sim¬ 
ulations, where instead of a constant current, a constant 
voltage drop is the control parameter and the current 
fluctuates around a saturated value [H, EE - 

From our experimental data, micro-vortices only set in 
from Regime III , slowly growing in size and speed with 
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FIG. 4: (a) The growth rate of vortex size (•) and voltage drop 
AV (►) during the linear convective Regime III indicated in 
Fig. 3, under different current densities, (b) The dependence 
of the average vortex size (Lmix) (B) and speed (0) on the 
applied current density across the membrane in the saturated 
convective regime IV. The error bar shows the standard devi¬ 
ation of the time averaged values, revealing more fluctuations 
in vortex dynamics at higher currents, (c) The average volt¬ 
age drops, AVoic, over the saturated vortex size in Regime 
IV, (•) obtained from the growth rate data of Regime III (a) 
and from the data of saturated AVji^ in Regime IV (•). (d) 
The average conductivity in the mixing layer in the growth 
regime, crm (◄), and in the saturated regime criv (•). 

time, accompanied by a linear increase of AV in time 
(e.g. 150 -450 s in Fig. [Hi). In addition, our experimental 
result of rms vortex velocity is quantitatively consistent 
with that found in 2D simulations of electro-osmotic in- 
stabilityunder similar electrical forcing (shown in Fig. 4c 
in Ref. [^), albeit different electrical boundary conditions 
(constant current vs. constant voltage). Based on these 
observations, the convective transport carried by swirling 
micro-vortexes is very likely the main cause of the over¬ 
limiting conductance observed. However, the fundamen¬ 
tal origin of the convective instability, which can be in¬ 
duced via an equilibrium or non-equilibrium mechanism 
as suggested by different theories [H, [^, [^, , remains 

elusive. Our data revealing quantitative growth of AVoic 
and vortex speed and size can motivate future theoreti¬ 
cal investigations under constant currents to identify the 
primary origin of the convective instability. 

We now show, for the first time, the dependence of 
electroconvective dynamics on the forced current density 
in Fig. m In Regime III, both the voltage AVoic and 
the mixing layer thickness of vortex region L^ix initially 
grow linearly (see Fig. |3^). These growth rates versus 


applied current densities are shown in Figure|4^. Fig.|4|o 
displays the average vortex speed and size in the satu¬ 
rated, over-limiting Regime IV, where unsteady dynam¬ 
ics are observed. Both vortex size and speed increase with 
current density, but the rms velocity increases stronger 
than the vortex size for the increasing current, underlin¬ 
ing the importance of convective transport in this regime. 
The large error bars at the higher current densities reflect 
the increasing fluctuations and chaotic motions of the in¬ 
dividual vortices. 

From the data of growth rates in Fig. [4^, the volt¬ 
age difference required to grow to a vortex region L^ix 
can be calculated using Regime III data: AV^/J = 
{6{AV)/St){6Lmix/6t)~^Lmix- Revealed in Fig. [H:, this 
voltage difference estimated for the final saturated vor¬ 
tex region agrees well with the voltage drop AVoic over 
the saturated mixing layers observed in Regime IV: 
AVji^ = AV — IRohm — AVc, indicating that micro¬ 
vortices only set in with an excess voltage drop, AVoic, 
across the mixing layer of vortex region. Furthermore, 
one could estimate the electrical conductivity a in this 
mixing region due to the presence of swirling vortices: 
tTiv = LmixJ/AVoi^, for each J. This conductivity is 
approximately constant for different current densities, 
and similar to the conductivity of the mixing layer in 
the growth regime III: crm = J{SLmix/St){S{AV)/St)~^, 
revealed in Fig. [Hi. We found that the conductivity 
in the mixing layer is similar for all the experiments: 
O’ = 55 ± 12/xS/cm. In contrast, in the under-limiting 
and limiting regimes, without vortices, ion concentra¬ 
tions are slowly depleted due to electrical migration and 
concentration diffusion, as well as the charge selectivity 
of the membrane. In the over-limiting regime, micro¬ 
vortices are initiated, continue to grow, and finally are 
saturated, right below the membrane. The viscous dissi¬ 
pation in the mixing layer, Pyisc, was estimated with the 
gradients of planar velocity field, using a nearest neigh¬ 
bor approach: Pyisc = M dV, with the Einstein 

notation {i,j = 1,2), volume element dV, and /r the liq¬ 
uid dynamic viscosity (1 mPa-s). We assume no shear 
in the 2 direction to estimate the volume integral of vis¬ 
cous dissipation. Comparing to the electrical power in¬ 
put Pin — Pelec ~ J.^mevaAVoic, the ratio of Pyisc/Pin id 
Regime IV is on the order of magnitude of 10“® — 10“®. 
This indicates extremely low power of energy dissipation 
by convective vortices, and thus AVoic represents the 
electrical resistance of the mixing layer of low ion con¬ 
centration. The advection flow modifies the anion con¬ 
centration gradient by bringing anions towards the mem¬ 
brane, which impairs the early-time ICP (in Regime II) 
and sustains over-limiting conductivity. 

In summary, the electrodialysis system involving a 
charge selective membrane for charge separation in gen¬ 
eral presents four distinct dynamic regimes observed in 
chronopotentiometric measurements: (I) a linear Ohmic 
electric response as charge diffusion and migration takes 
place; (II) a jump in electrical response (AV) during the 
development of ICP due to charge-selectivity of the mem- 
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brane; (III) a linear growth regime where micro-vortices 
grow in both size and speed with time; and finally (IV) a 
saturated electro-convective regime having saturated val¬ 
ues of vortex speed/size and voltage response AV. Our 
quantitative results of the growth rates and saturated 
electro-convective responses elucidate that micro-vortices 
only set in with an excess voltage AVoic, have small vis¬ 
cous dissipation, and moreover sustain a nearly constant 
conductivity in the mixing region. 
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